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We study a Centauro-type phenomenon in high-energy heavy-ion colhsions by 

Q^ ■ assuming that pions are produced semiclassically both directly and in pairs through 

(N' 

m 



the isovector channel 



The leading-particle effect and the factorization property of the scattering 

'NT ■ 

^ I amplitude in the impact-parameter space are used to define the classical pion field. 

p ■ By analyzing the joint probability function Pii^{nQ,n_) for producing uq neutral 

CLi' and n_ negative pions from a definite isospin state JJ3 of the incoming leading- 

^ ■ particle system we show that only direct production of pions without isovector pairs 



favors Centauro-type behavior. The presence of isovector pairs seems to destroy the 
effect. Our conclusion is supported through the calculation of two pion correlation 
parameters, /2~ and /2°, and the average number of neutral pions ((?^o)n_) as a 
function of negative pions (n_) produced. 
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In recent years several cosmic-ray experiments [1] have reported evidence for 
the existence of Centauro events characterized by an anomalously large number 
of charged pions in comparison with the number of neutral pions, indicating that 
there should exist a very strong long-range correlation between two types of the 
pions. Such long-range correlations are possible if pions are produced coherently 
and constrained by the global conservation of isospin [2-7] . 

Although the actual dynamical mechanism of the production of a classical pion 
field in the course of a high-energy collisions is not known, there exist a number 
of interesting theoretical speculations [8-13] that localized regions of misaligned 
vacuum might occur in ultrahigh-energy hadronic and heavy-ion collisions. These 
regions become coherent sources of a classical pion field. In early models [2,3], 
however, the coherent production of pions is taken for granted and is considered to 
be a dominant mechanism. 

These models also predict strong negative correlations between the number 
of neutral and charged pions. In fact, the exact conservation of isospin in a pion 
uncorrelated jet model is known [4,5] to give the same pattern of neutral/charged 
fluctuations as observed in Centauro events. This strong negative neutral-charged 
correlation is believed to be a general property of the direct pion emission in which 
the cluster formation ( or the short-range correlation between pions ) is not taken 
into account [6,7,14]. 

A similar conclusion may be drawn when a proper multipion symmetrization 
combined with the isospin singlet channel is considered [15]. 

In this paper we consider the leading-particle effect as a possible source of a 
classical pion field [16]. Pions are assumed to be produced from a definite isospin 
state of the incoming leading-particle system both directly and in isovector pairs. 



Coherently emitted isovector clusters decay subsequently into pions outside the 
region of interaction. We discuss the limiting behavior of the probability distribution 
of neutral pions, the neutral - to - charged ratio and the corresponding two-pion 
correlation function /2^, as well as the variation of the average number of neutral 
pions ((no)n_) as a function of the number of negative pions (n_) produced. 

We show that the Centauro effect is strongly suppressed if pions are produced 
in isovector pairs. At high energies most of the pions are produced in the cen- 
tral region. To isolate the central production, we adopt high-energy longitudinally 
dominated kinematics, with two leading particles retaining a large fraction of their 
incident momenta. The energy available for the hadron production is 

Ehad = -^yS — Eieading- (1) 

It is related to the inelasticity K of the collisions by 

Ehad = ^V^K. (2) 

Note that at fixed total cm. energy ^/s, the hadronic energy Ehad varies from event 
to event. 

With a set of independent variables s, {qir, Hi} = Qi, i = 1,2, . . .n, the n-pion 
contribution to the s-channel unitarity becomes an integral over the relative impact 
parameter b of the two incident leading particles: 

1 /■ " 
Fn{s) = -r d^bll dqi I T„(s, h-l...n)\\ (3) 

where dq = d'^grdy/ 2 (2,71)^. The normalization is such that 

F„(s) = sarXs), 

00 

(Tineli-S) = J^'^nis)- (4) 

n=l 



At high energies the matter distributions of the two-leading-particle system are 
Lorentz-contracted disks in the cm. system. We assume that, in the early stages 
of the collision, a highly excited localized system occurs that relaxes through the 
coherent emission of pions [17]. According to [18], this type of coherent emission of 
pions should saturate close to the threshold since, once it starts, the resulting pion 
emission, in the absence of a " resonance cavity " , prevents further buildup of pion 
fields. The basic assumption of the independent pion-emission model, neglecting the 
isospin for a moment, is the factorization of the scattering amplitude Tn{s, 6; 1 ... n) 
in the b space: 

^ -n-l n 

T„(s, b;l...n) = 2sfis, b)^=l[ J(s, b; g,), (5) 

where, owing to unitarity, 

|/(s,6)p=e-^(^'^^), (6) 

and 

n{s,b)=ldq\J{s,b;q)\' (7) 

denotes the average number of emitted pions at a given impact parameter b. The 
function | J{s,b;q) p, after the integration over b, controls the shape of the single- 
particle inclusive distribution. A suitable choice of this function also guarantees that 
the energy and the momentum are conserved on the average during the collision. 

The above particular factorization of T„ was considered earlier by R. Avin 
et al. [19] in connection with the construction of a unitary model of multiparticle 
production. 

The inclusion of isospin in this model is straightforward [3] . We observe that 
the factorization of T„ is a consequence of the pion field satisfying the equation of 
motion 

{n + fj,^)7f{s,b;x)=j{s,b;x), (8) 



where j is a classical source related to J{s, b; q) via the Fourier transform 

J{s, b;q)= f d^xe^'^^'jis, b; x). (9) 



The standard solution of Eq.(8) is given in terms of in- and out-fields that are 
connected by the unitary S matrix S{b, s) as follows: 

where 

TTclassical = / d^x' A{x - X] fl)j{s, b] x'). (11) 



The S matrix following from such a classical source is still an operator in the 
space of pions. Inclusion of isospin requires S{s, b) to be also a matrix in the isospace 
of the leading particles. 

The coherent production of isovector clusters of pions is described by the fol- 
lowing S matrix: 

Sis,b)= fd^e\e)DiY..l;s,b){e\, (12) 

where | e ) represents the isospin-state vector of the two-leading-particle system. The 
quantity D{J2cJc',s,b) is the unitary coherent-state displacement operator defined 
as 

D{Y^ X; s, b) = explY^ / dqJc{s, b; q)d'J{q) - H.c], (13) 

c c '' 

where ac (q) is the creation operator of a cluster c and the summation J2c is over all 
clusters. The clusters decay independently into c = 1,2,... pions outside the region 
of strong interactions. 

Clusters decaying into two or more pions simulate a short-range correlation 
between pions and behave as p, p', ■ ■ ■, etc. in isospace. They need not be well- 
defined resonances. The more pions in a cluster, the larger the correlation effect 



expected. If the conservation of isospin is a global property of the colliding system, 
restricted only by the relation 

/ = /'+/;, (14) 

where I^, denotes the isospin of the emitted pion cloud, then Jc{s, b; q) is of the form 

Jc{s, b; q) = Jc{s, b; q)e, (15) 

where e is a fixed unit vector in isospace independent of q. The global conservation 
of isospin thus introduces the long-range correlation between the emitted pions. 

If the isospin of the system of two incoming ( outgoing ) leading particles is 1/3 
and J'/g, respectively, then the initial-state vector of the pion field is S{s,b) \ II3), 
where | II3) is a vacuum state with no pions but with two leading particles in the 
isostate characterized by II3. Then the n-pion production amplitude is 

iTn{s, b;qi... g„) = 2s{I'l'^] gi . . . g„ | S{s, b) \ Ih). (16) 

The unnormalized probability distribution of producing n+vr^, n_7r-, and n^Ti^ pions 
is defined as 

W{n+nnQ,I'I'^Jh) = j d%dqidq2 . . . dqn \ {I'l'^n+njiQ \ S{s,b) \ J/g) p, (17) 
where n = n+ + n_ + no- 

Assuming further that all (/', Jg) are produced with equal probability, we can 
sum over all possible isospin states of the outgoing leading particles to obtain 



^ , , ErnW{n+n_no,rr,;Ih) 

P//3(n+n_no) 



E„+„_„o T.i'1^ W{n+n_no, /'/^; II3) 
This is our basic relation for calculating various pion-multiplicity distributions, pion 



multiplicities, and pion correlations between definite charge combinations. ^ 

In order to obtain some more detailed results for multiplicity distributions and 
correlations, one should have an explicit form for the source function Jc{s, b]q),c = 
1,2, 



'-! ■"! 



The results on the isospin structure are most easily analyzed in the so-called 
grey- disk model in which 

n,{s,b)=n,is)0{bois)-b), (19) 

where nc{s) denotes the mean number of clusters of the type c, bQ{s) is related to 
the total inelastic cross section, and 6' is a step function. 

Let us assume that pions are produced both directly and through isovector 
clusters of the p type. In this case, the distribution of neutral pions becomes 

2 (/ + /3)! J-1 no! 

where A{x) = n^x^ + (1 — x'^)np. Here n^ denotes the average number of directly 
produced pions, and rip denotes the average number of p-type clusters which decay 
into two short-range correlated pions. The function Pj^{x) denotes the associate 
Legendre polynomial. 

Note that the total number of emitted pions is 

(n) =n^ + 2np. (21) 

The correlation /2 between the pions, 

h = {n{n - 1)) - {nf = 2n„ (22) 



^In general, the probability W{n+n_nQ,I' l!^]II^) depends on {I'l'^^) dynamically. The final- 
leading-particle production mechanism usually tends to favor the (/', /a) ~ (/, /s) case, for example, 
if the final-leading particles are nucleons or isobars. However, if the leading particles are colliding 
nuclei, it is reasonable to assume almost equal probability for various {I\ I^) owing to the large 
number of possible leading isobars in the final state. 



is positive indicating deviation from the Poisson distribution (/2 = 0) owing to the 
emission of p-type clusters. 

In Fig. 1 we show the behavior of Pu^ino) for (n) = 50 and different combina- 
tions of (n^, rip) when (a) / = J3 = 0, and ( b) / = /a = 1. We see that Centauro-type 
behavior is obtained only for pions that are produced directly, i.e., for n^r 7^ and 
Up = 0. 

If i? = jA- denotes the fraction of neutral pions, then it is easy to see that in 
the limit (n) -^ 00 with R fixed, the probability distribution {n)Pjj,^{nQ) scales to 
the limiting behavior: 

{n)Pn,{no) - Pn,{R) = (/ + ^7^— ^ I Pi'i^n) \" [(1 " 37p)(i? - IpTK (23) 



2'{I + l3)\ 



where Ip = j;^ 



Up 

in)- 


= 7= — , ^0— \ aiid 






^2 _ R-ip 

l-37p 


0<7,<i 



(24) 

This limiting probability distribution is different from the usual Gaussian random 
distribution for which one expects peaking at i? = | as {n) -^ 00. 

Instead, we find peaking at i? = 7^ although 

(«o)// = {n) - {nch)ii 
1 



2/ + 3 



n^ + 2(/ + l)np] (25) 



and (no) = -{nch) = -(^); for either / = or ^tt = Up. 

If 7p = |, that is if n.,^ = fip, the probability Pii^{no) becomes a pure Poisson 
distribution and is independent of II3. 



In Fig. 2 we show the average number of neutral pions ((^o)n_)) as a function 
of the number of negative pions (n_) produced for different pairs of (n,r,^p) and 
I = I3 = 1. We note that {no)n_ decreases with n_ only in the region where nj, ^ Tip 
and (n) ^ n_, indicating again that the Centauro-type effect cannot be expected if a 
significant number of p-type clusters of pions is produced. Similarly we find that the 
neutral-pion two-particle correlation parameters /2" = {noUa) — (^o)(^a) — (^o)^oa 
for a = and — are related as 

ftn + lf2'n = n,^, (26) 



where 



4(J + 1) 



^2.^^" (2/ + 3)2(2/ + 5)^'''^ ""'^ ' ^^^^ 

showing that /^ 7/ is negative if rip is small. We should also mention earlier works 
on neutral-pion correlation effects [20]. 

In summary, the results of the present analysis show that the emission of 
isovector clusters of pions in the framework of a unitary eikonal model with the 
global conservation of isospin strongly suppreses the Centauro-type effect for pions. 
This might suggest that the Centauro-type effect, if any, will probably appear only 
in very limited regions of phase space where isovector clusters should be missing. 
How it is possible dynamically is not yet clear [10-13]. Similar discussion about 
emission of isoscalar clusters of pions is presented in [15]. 
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Figure captions : 

Fig. 1. Multiplicity distributions Pii,^{no) of neutral pions for (n) = 50 when 
the total isospin of the incoming-leading-particle system is (a) / = Is = and 
(6) / = Is = 1. The curve s represent different combinations of (nT^jTip), the aver- 
age number of directly produced pions, and the average number of p-type clusters, 
respectively. 

Fig. 2. The average number of neutral pions as a function of the number of 
negative pions produced for I = I^ = 1. The curves represent different combinations 
of {n^,np). 
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